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摘 要 采用 FLUENT 中 8 个 满 流 模型 以 及 曲率 修正 后 的 模型 对 U 型 弯 管 中 的 二 维 流动 进行 了 模拟 ， 详 细 对 比 
了 两 个 雷诺 数 105 和 105 下 速度 、 涪 动能 、 湛 流 剪 切 应 力 、 摩 擦 系数 和 静 压力 系数 的 实验 和 模拟 结果 .研究 发 现 ， 当 
Re = 105，Realizable k 一 模型 、Spalart-Allmaras 模型 、RNG k 一 。 模型 能 准确 模拟 出 流动 分 离 的 位 置 ; 当 Re=106， 
所 有 模型 预测 到 的 分 离 位 置 都 滞后 . 曲率 修正 对 SKE 以 外 的 灌流 模型 预测 性 能 基本 上 没有 改善 ;SKE 模型 对 曲率 修正 的 
影响 最 为 敏感 , 对 速度 、 满 动能 、 内 壁面 处 的 摩擦 系数 的 模拟 结果 都 有 一 定 的 改善 . 
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Abstract The flow in 2D of Sha is simulated, with eight common turbulence models and 
curvature correction model oe By comparing the simulation and experimental results of 
velocity, turbulence kin ergy, turbulent shear stress, friction coefficient and static pressure 
coefficient under two kiids’of Reynolds number 105 and 10°, it is find that the Realizable 大 一 s model, 
Spalart-Allmaras model and RNG k 一 = model have a good simulatiCpsresults about the position of 
the flow separation when the Reynolds number is 105. BesidegS he separation positions predicted 
by all models are lagging behind. The models with curvaturé\correction are not better than before 
expect SKE. The model of SKE is more sensitive to the Cnce of curvature, and there is a greater 
improvement on simulation about velocity and turb ene kinetic energy:. 
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0 引言 < 


航空 发 动机 和 燃气 轮机 (简称 “两 机 ”) 应 用 广 。 类 是 雷诺 应 力 模型 (RSM). 尽管 RSM 模型 能 够 模 
泛 , 发 展 先进 “两 机 ” 技术 , 对 动力 、 节 能 、 减 排 、 降 拟 出 一 些 局 部 典型 复杂 流动 现象 , 但 是 由 于 仍然 留 
噪 等 都 有 重大 意义 . 设计 高 性 能 “两 机 ”"， 离 不 开 先 。 有 固有 的 物理 不 足 ， 以 及 收敛 性 不 好 , 较 大 的 复杂 
进 的 CFD 手段 、 目 前 ， 准 确 预测 “两 机 ”内 的 复杂 ”性 和 计算 量 等 多 种 原因 ， 至 今 难 以 应 用 于 航空 发 动 
流动 (如 大 太 度 分 离 等 ) 对 当前 CFD 提出 了 巨大 的 。 机 的 研制 . 而 涡 务 滑 流 模型 都 是 基于 人 简单 边界 、 平 衡 
挑战 . DNS、LES、hybrid LES/RANS 等 高 精度 方法 汕 流 、 各 向 同性 假设 发 展 起 来 的 . 航空 发 动机 中 的 流 
由 于 计算 耗费 太 大 而 无 法 应 用 于 日 常 工程 设计 , 而 动 非常 复杂 , 存在 着 大 分 离 流 动 和 复杂 涡 系 , 常 处 于 
工程 广泛 应 用 的 RANS 方法 存在 灌流 模型 模拟 准确 。 洪流 非 平 衡 、 各 向 异性 的 状态 , 使 得 常用 涡 秋 模型 难 
性 的 间 题 中 以 准确 模拟 . 大 曲率 是 航空 发 动机 中 流动 复杂 的 主 

现 有 汗 流 模型 主要 分 为 两 大 类 : 一 类 是 以 ”要 影响 因素 之 一 ， We madi 
Boussinesq 涡 务 性 假设 为 基础 的 涡 秋 性 模型 ， 另 一 “” 率 对 流动 的 影响 . 在 这 种 背景 下 ,Spalart 和 ShurD 
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提出 了 一 种 曲率 修正 方法 ， 并 有 效 提 高 了 Spalart- 
Allmaras 模型 的 精度 . Smirnov 和 MenterBl 把 这 一 
修正 方法 应 用 到 了 SST k-w 模型 中 ， 也 有 效 提高 了 
该 模型 的 精度 .但 到 目前 为 止 ， 对 大 曲率 流动 中 各 
用 清流 模型 性 能 的 系统 研究 还 较 少 ， 亚 需 进一步 开 
展 深 入 研究 。 

U 型 弯 管 流动 很 好 地 体现 了 大 曲率 对 流动 的 影 
啊 ， 因 此 成 为 了 研究 曲率 影响 的 经 典 流 动 。 其 中 
NASA 的 Monson 等 85 对 U 型 弯 管 流动 进行 
了 详细 的 实验 测量 ， 提 供 了 丰富 的 流 场 数据 。 本 文 
针对 该 二 维 U 型 弯 管 流动 ， 采 用 最 常用 的 8 个 
庙 流 模型 对 其 进行 数值 模拟 . 8 种 清流 模型 分 别 是 
Spalart-Allmaras 模型 [6 (SA)、Standard 有 一。 模 
型 门 (SKE)、Realizable 大 -= 模型 [S| (RKE)、RNG 
k 一 es 模型 四 (RNG KE)、Standard k-w 模型 60| 
(SKW)、SST 一 ww 模型 0113(SST KW)、v2-f 模 
型 5544(V2F) 和 雷诺 应 力 模型 05,19(RSM)。 通过 
模拟 结果 与 实验 结果 的 对 比 ， 系 统 的 分 析 各 济 流 模 
型 对 大 曲率 流动 的 预测 能 力 。 进 一 步 采 用 带 有 
修正 的 SA、SKE. RNG KE、RKE、SKW 和 8 
对 该 算 例 进行 了 模拟 ， 分 析 曲 率 修 正 对 误 
全 hj 辐 [I 
能 的 影响 。 EN 
1 二 维 U 型 弯 管 模型 
1.1 二 维 型 栓 管 参数 向 站 

本 文 模拟 的 二 维 U 型 弯 管 参数 以 及 实验 数据 均 
来 自 Monson 等 bh17 的 实验 , U 型 管 参数 如 图 1 所 
示 . 管 高 度 五 = 3.81 cm， 入 口 段 长 度 为 10 万 = 38.1 
cm， 出 口 段 长 度 为 12 万 = 45.72 cm， 弯 管 处 内 壁 半 
径 mr= 1.91 cm, 外壁 半 径 ro= 5.72 cm, s 为 沿 U 型 


管 中 心 线 的 自然 坐标 系 ， 以 入 口 段 和 弯 管 交界 处 以 J" 


原点 , 入 口 段 为 负 , 弯 管 处 和 出 口 段 为 正 . 了 辆 欧 蔷 
标 值 在 内 壁面 处 为 零 ， 且 恒 为 正 值 ， 外 壁面 处 最 大 
值 为 互 , X 轴 坐 标 值 在 入 口 段 为 负 , 出 口 段 为 正 , 弯 
管 处 采用 极 坐 标 ， 原 点 为 内 外 壁面 的 圆心 。 作 图 时 
位 置 都 用 管 径 互 做 了 无 量 纲 化 处 理 , 速度 都 以 入 口 
速度 Uior 为 基准 做 了 无 量 纲 化 处 理 ， 
1.2 实验 数据 

实验 数据 来 源 于 Monson 等 向 的 实验 , 进口 速 
度 为 et= 32 m/s (Ma 0.1)， 进 口 总 压 分 别 为 
P=1.2Po 和 12 记 , 忆 为 标准 大 气压 ,出口 压 力 与 
进口 压力 一 样 , 分 别 可 以 达到 进口 雷诺 为 Re = 105 
和 105, 实验 数据 测量 主要 运用 二 维 LDV (激光 多 普 
勒 测速 ) 技术 测量 了 流动 速度 ， 此 外 还 测量 了 U 型 
管 的 表面 摩擦 力 、 汗 动能、 涡流 切 应 力 压力 等 参数 ， 
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关于 更 多 的 数据 测量 细节 可 以 参考 Monson 等 向 的 
工作 . 
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图 1 U 型 弯 管 参数 示意 图 外 
Fig. 1 The geometry of U-bend 国 


2 类 算 设 置 
从 网 阁 的 设置 

XT 为 了 减 小 网 格 对 模拟 结果 的 影响 ,根据 各 滑 流 
模型 对 壁面 网 格 的 要 求 ， 在 近 壁面 采用 的 增强 壁面 
模型 , 近 壁 面 网 格 y+ ~1. 同时 本 文通 过 调整 网 格 分 
布 以 及 网 格 数 ,采用 了 三 套 网 格 , 如 图 2 和 表 1 所 
示 , 进行 了 网 格 无 关 性 检验 . 通过 模拟 进口 Re = 105 
条 件 下 不 同 截面 的 速度 分 布 曲线 来 对 网 格 无 关 性 进 
行 分 析 . 


图 2 U 型 弯 管 网 格 示 示 意 外 
Fig. 2 The mesh of U-bend 


表 1 网 格 节点 分 布 情况 
Table 1 The grid node distribution 


节点 数 a b c h 网 格 总 数 
第 一 套 104 125 125 97 34388 
第 二 套 125 150 150 117 49735 
第 三 套 150 180 180 140 71400 


三 套 网 格 均 为 SKE 模型 对 进口 Re = 105 的 
二 维 U 型 弯 管 进行 模拟 .通过 对 s/H = -4、0 = 
0°、9 一 90°、9 二 180°、s/ 瑟 二 2、s/ 且 ==4 以 及 出 口 
位 置 的 速度 分 布 进行 对 比 (限于 篇 幅 没 有 给 出 )， 发 
现 三 套 网 格 的 模拟 结果 完全 吻合 ， 做 到 了 网 格 无 关 
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解 。 因 此 ， 本 文采 用 第 二 套 网 格 进行 后 期 的 数值 模 
拟 计算 .。 
2.2 边界 条 件 及 满 流 模 型 

根据 实验 文献 由 中 提供 的 数据 , 选取 上 游 处 距 
离 弯 曲 段 进口 45 处 的 物理 量 作 为 进口 条 件 。 马赫 
数 为 0.1 时 静 压 与 总 压 的 比值 x = 0.993. 雷诺 数 Re 
= 105 时 , 进口 静 压 设置 为 P = 120787 Pa; 当 雷 诺 
数 Re = 105 时 , 进口 静 压 设置 为 已 = 1207870 Pa. 
出 口 设置 为 压力 出 口 静 压 与 进口 静 压 相等 。 

为 了 比较 不 同 灌流 模型 对 模拟 结果 的 影响 ， 在 
本 文 的 研究 中 , 采用 SA 模型 、SKE 模型 、RKE 模 
型 、 RNG KE 模型 、 SKW 模型 、SSK KW 模型 、V2F 
模型 和 RSM 模型 8 种 常用 的 清流 模型 对 U 型 管内 
的 流动 进行 了 计算 ，SKE 模型 、RKE 模型 、RNG 
KE 模型 和 RSM 模型 在 近 壁 区 采用 了 增强 壁面 处 理 
模型 . 为 了 研究 曲率 修正 后 的 灌流 模型 的 模拟 性 能 ， 
本 文 最 后 也 对 带 曲 率 修正 的 SA CC 模型 、SKE CC 
模型 、 RKE CC 模型 、 RNG KE CC 模型 、 SKW CC 
模型 和 SST KW CC 模型 对 流动 进行 了 模拟 . 


夏 


3 常用 济 流 模型 的 模拟 结 
3.1 速度 模拟 结果 

从 图 3 (a) 和 图 3 (b) 中 可 以 看 到 ,由 于 曲率 影 
响 ，0? 截面 位 置 流 速 已 经 发 生 了 偏 移 ， 由 于 气流 
流动 中 产生 的 离心 力 与 径 向 压力 梯度 的 不 平衡 ， 导 
致 U 型 管 中 形 成 一 个 较 大 的 低压 区 ， 从 而 使 其 流 
速 增加 ， 而 外 侧 壁 面 由 于 较 大 离心 力 的 作用 ， 使 
外 侧 壁 面 处 产生 一 个 较 大 的 高 压 区 ， 使 得 外 侧 璧 
面 处 速度 较 低 。 在 该 截面 处 ， 所 有 模型 在 内 侧 璧 
面 处 的 模拟 结果 都 偏 大 ， 各 油 流 模型 模拟 结果 相差 
不 大 . 

在 90° 截面 位 置 ， 内 侧 速 度 达到 了 最 大 值 ， 外 
侧 速度 降 为 最 低 值 ， 如 图 3 (c) 和 图 3 (d)。 该 截面 
处 ,内 外 人 币 殿 面 的 模拟 结果 同 实验 结果 有 明显 的 差 
异 。 SM 模型 在 外 侧 壁 面 处 的 模拟 结果 模拟 

，SKW 模型 和 SST KW 模型 在 内 侧 壁 面 处 

拟 结果 同 实验 吻合 的 最 好 . 整体 来 说 ， RSM 模 


:OD 型 模拟 的 结果 同 实验 最 为 接近 . 
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图 3 弯曲 段 各 截面 的 速度 分 布 曲线 


Fig. 3 The cross section of velocity distribution 


在 180° 截面 位 置 , 各 满 流 模型 模拟 结果 差异 较 和 出 所 


大 , 如 图 3 (e) 和 图 3 (f). 当 Re = 105, RKE 模型 的 
模拟 结果 与 实验 结果 最 接近 , 而 当 Re = 105, RSM 
模型 与 实验 结果 最 为 接近 。 


SN\ 
从 图 3 可 以 看 到 , 相 比 Re = 105， 各 满 } 
在 Re = 10 时 的 模拟 结果 与 实验 结果 更 台 > 


We 
性 , 本 文 对 Re = 10° 流动 下 , 90* 瑚 次 位 置 的 消 动 能 
的 分 布 情况 进行 了 进一步 研究 , 铬 狠 4 所 示 , 该 截面 
在 外 侧 壁 面 附近 滑动 能 的 灾 验 绪 果 出 现 了 峰值 ,各 
清流 模型 的 模拟 结果 偏 小 . 壬 合 图 3 (c) 和 图 3 (d) 
中 该 截面 的 速度 分 布 情况 ， 外 侧 壁面 处 速度 模拟 结 
果 偏 小 可 能 是 由 于 该 处 滑动 能 模拟 结果 不 准确 造成 
的 . 在 内 侧 壁 面 处 , SKW 和 SST KW 模型 出 现 了 峰 
值 , 这 与 实验 结果 偏差 较 大 . 整体 来 看 , 对 于 滑动 能 
的 模拟 结果 , RSM 模型 模拟 结果 更 加 吻合 . 

图 5 显示 了 雷诺 数 Re = 105 和 Re 二 10 于 
下 游 5/ 了 = 2、s/H= 4 和 出 口 截 面 处 流速 分 贸 胸 模 
拟 情况 . 图 5 显示 , 相对 于 弯曲 段 , 在 下 游 直 管 段 不 
同市 流 模型 问 的 模拟 结果 差异 性 较 大 ， 各 消 流 模型 
模拟 结果 与 实验 结果 相差 较 大 . 从 图 5 中 可 以 看 到 ， 
在 内 侧 壁 面 处 各 模型 模拟 结果 都 偏 小 ， 随 着 流动 向 
下 游 发 展 ,这 种 偏差 在 逐渐 减 小 , 如 图 5 (e) 和 图 5 
(f)， 在 出 口内 侧 壁 面 处 模拟 结果 与 实验 结果 的 偏差 
较 小 . 从 图 5 (c) 和 图 5 (d) 可 以 看 到 ,与 实验 结果 
相 比 ， 模 拟 结果 得 到 了 一 个 较 薄 的 边界 层 ， 而 在 出 
口 处 , 如 图 5 (e) 和 图 5 (9， 当 流动 逐渐 从 曲率 影响 
中 恢复 时 ， 滑 流 模型 预测 的 结果 恢复 的 较 慢 . 

为 了 进一步 分 析 下 游 的 直 管 流动 情况 ,图 6 和 
图 7 中 给 出 了 雷诺 数 Re = 106 时 , 在 s/ 矿 =?2 


立 置 的 清 动 能 模拟 结果 .从 图 6 中 可 
在 s/ 瑟 = 2 截面 ， 内 侧 壁 面 附近 模拟 的 
能 峰值 均 比 实验 值 偏 小 ， 而 在 弯 管 中 心 靠近 外 
则 壁面 附近 出 现 了 应 动能 平台 区 ， 与 实验 结果 有 所 
偏差 。 


以 


30 40 50 60 70 
103K/ U2 
图 4 湛 动 能 分 布 图 9 = 90°, Re = 105 
Fig. 4 Turbulent kinetic energy distribution 
0 = 90°, Re = 10° 


结合 图 5 (a) 和 图 5 (b), 在 相同 位 置 出 现 了 速度 
平台 以 及 在 内 侧 壁 面 处 模拟 结果 偏 小 ， 可 能 是 由 于 
各 水 流 模型 对 涡 动 能 的 模拟 结果 不 准确 导致 的 . 其 
中 SKE 模型 的 模拟 结果 偏差 最 大 ,而 RSM 模型 模 
拟 结果 吻合 较 好 。 从 图 7 中 出 口 消 动 能 实验 结果 显 
示 ， 沸 动能 已 逐渐 恢复 平稳 ， 而 各 模型 的 济 动 能 模 
拟 结果 仍然 有 较 大 的 波动 ,因此 从 另 一 个 角度 也 说 
明 满 流 模型 预测 的 流动 从 曲率 影响 中 恢复 比 实验 结 
果 更 晚 一 些 . 
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图 5 下 游 直 管 付 各 帘 面 速度 分 布 图 


Fig. 5 The velocitRgia s of the downstream sections 
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图 6 s/H =2 处 湾 动 能 分 布 Re = 105 图 7 出 口 处 满 动能 分 布 Re = 105 
Fig. 6 Turbulent kinetic energy distribution s/H = 2, Fig. 7 Turbulent kinetic energy distribution Outlet, 
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3.2 分 离 区 的 模拟 结果 

根据 Sandbornl' 引 的 弯 管 实验 可 以 发 现 ， 在 雷 
诺 数 Re = 105 时 ,分离 位 置 在 170° 左右 , 分 离 区 
平均 最 大 高 度 为 0.04 互 , 分 离 再 附 的 位 置 在 下 游 0.5 
H~18 的 位 置 处 . 在 雷诺 数 Re = 105 时 , 通过 实验 
发 现 , 分 离 位 置 在 150° 位 置 , 分 离 区 平均 最 大 高 度 
为 0.14 五 , 分 离 再 附 的 位 置 在 下 游 1 五 ~~1.5 玉 的 
位 置 处 。 

表 2 和 表 3 显示 了 两 种 流动 条 件 下 ， 各 泣 流 模 
型 对 流动 分 离 的 模拟 情况 


表 2 Re = 10” 时, 分 离 位 置 统 计 情 况 


Table 2 The separation position statistics, 


Re = 105 

模型 /实验 ”分离 位 置 /(? ) 再 附 位 置 竖 直 高 度 

实验 170 0.5H ~1.5H 0.04H ~ 0.16H 

SA 160 1.67H 0.016H 

SKE 无 无 无 

RKE 160 1.23H 0.11H 
RNG KE 162 1.15H 0.08H 

SKW 145 1.71H 0.19H 如 
SST KW 140 1.98H 0.24 A 

V2F 155 2 万 


RSM 128 


表 3 Re= 10 时 , 分 离 计 情 况 
Table 3 The 0 ition statistics, 


Re = 

模型 /实验 ”分 离 位 置 /(°) 再 附 位 置 竖 直 高 度 

实验 150 108 ~15H8 0148 ~ 0.238 

SA 171 1.4H 0.09H 

SKE 无 无 无 

RKE 160 1.36H 0.14H 
RNG KE 165 1.26H 0.12H 

SKW 164 1.2H 0.13H , 
SST KW 162 1.47H 0.16H x 

V2F 166 1.52H 0.17 

RSM 162 1.17H 0.16 


从 表 2 中 , 可 以 看 到 当 Re = 105 时 , SKE 模型 模 
拟 结果 是 最 差 的 , 没有 模拟 到 流动 的 分 离 情 况 . RKE 
模型 和 RNG KE 模型 对 分 离 点 位 置 模拟 最 接近 实验 
分 离 点 , 而 RSM 模型 的 分 离 位 置 要 明显 比 实验 结果 
提前 . 而 对 于 再 附 位 置 ， SKW 模型 、SST KW 模型 
和 V2F 模型 比 实验 结果 推迟 了 很 多 , 这 些 模型 对 分 
离 涡 高 度 的 模拟 结果 也 偏 大 。 整体 来 看 ,SKW 模型 
和 SST KW 模拟 到 的 分 离 区 更 大 一 些 . 各 模型 模拟 
分 离 位 置 点 都 比 实验 结果 靠 前 ,再 附 位 置 比 实验 结 
果 推 后 ,分 离 区 平均 最 大 高 度 也 较 大 ， 因 此 对 于 分 
离 涡 的 大 小 ,整体 模拟 结果 比 实验 结果 要 大 . 

当 Re = 106 时 , 如 表 3 所 示 ,， 各 模型 对 分 离 位 


置 的 模拟 结果 都 比 实验 结果 推 后 ， 其 中 SA 模型 偏 
差 最 大 . 对 比 表 2 和 表 3, 各 消 流 模型 在 较 高 的 雷诺 
数 时 , 对 分 离 涡 的 模拟 更 加 准确 , 而 SKE 模型 在 两 
种 流动 条 件 下 ， 模 拟 结果 均 最 差 ， 没 有 模拟 到 分 离 
流动 的 存在 . 整体 来 看 ， 当 雷诺 数 较 高 时 ，RSM 模 
型 和 RNG KE 模型 对 分 离 的 预测 能 力 较 好 . 
3.3 表面 摩擦 系数 和 壁面 静 压力 系数 的 模拟 结果 
图 8 和 图 9 中 显示 了 两 种 流动 情况 下 ， 内 外 侧 
壁面 的 摩擦 系数 绝对 值 曲线 。 


©® ©®EXP Re=105 


(b) 外 壁面 
图 8 内 外 壁面 的 摩擦 系数 绝对 值 分 布 Re = 105 


Fig. 8 The wall friction coefficient absolute value distribution 
of inner wall and outer wall, Re = 105 


当 流 动 接近 弯曲 段 时 , 内 侧 壁 面 流 体 加 速 , 摩擦 
系数 Cr 急剧 增加 . 在 Re = 10° 时, 各 消 流 模型 模拟 
到 弯曲 段 的 内 侧 壁 面 的 摩擦 系数 峰值 明显 偏 大 ， 如 
图 8 (a). 而 在 Re = 10" 时 , 该 处 的 模拟 结果 有 所 改 
善 , 如 图 9 (a). 在 弯曲 段 后 部 ， 流 动 减速 ，Cf 快速 
降低 ,发 生 分 离 , 出 现 了 零 值 点 .从 图 8 (a) 中 也 可 
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以 看 到 , 在 Re = 105 时 , 各 潮流 模型 预测 的 再 附 点 
都 比 实验 结果 推 后 ， 导 致 下 游 摩擦 系数 偏 小 。 而 在 
Re = 106 时 , 各 模拟 结果 模拟 到 的 再 附 点 同 实验 结 
果 吻 合 较 好 ， 如 图 9 (a)。 


(b) 外 壁面 


图 9 内 外 壁面 的 摩擦 系数 绝对 值 分 布 Re = 105 QQ 


Fig. 9 The wall friction coefficient absolute value dist ia 
of inner wall and outer wall, Re = 106 NS 


除了 摩擦 系数 Ct 以 外 ,壁面 静 压 力 系 数 Cb 在 
工程 中 也 是 非常 重要 的 参数 . 图 10 和 图 11 显示 
了 各 灌流 模型 的 壁面 更 压 力 系数 Cb 绝对 值 的 模拟 
结果 。 

同 摩 擦 系数 相反 的 ， 如 图 10 (a) 和 图 11(a) 所 
示 , Cp 在 弯曲 段 内 侧 壁 面 , 流速 增加 , 压力 降低 , 更 
压 系数 达到 最 低 的 峰值 ， 随 着 流速 的 下 降 ， 压 力 降 
到 最 低 值 后 ,急剧 上 升 , 产生 很 大 的 径 向 压力 梯度 ， 
导致 流动 分 离 ， 最 终 曲线 逐渐 恢复 至 平 线 。 由 于 在 
此 过 程 中 ， 有 流动 损失 形成 ， 因 此 下 游 的 静 压 力 系 
数值 要 低 于 进口 段 的 值 . 同时 ， 从 图 10 以 看 到 , 在 
分 离 区 之 前 ， 各 注 流 模型 对 压力 系数 的 模拟 结果 差 
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异性 很 小 ， 而 且 同 实验 吻合 程度 较 大 ， 而 当 流动 分 
离 后 ， 各 清流 型 对 流动 模拟 有 明显 的 不 同 . 在 下 游 
分 离 后 ， 各 消 流 模型 模拟 结果 都 比 实验 偏 大 . 整体 
来 看 ， RSM 模拟 结果 相对 较 好 。 


(b) 外 壁面 


图 10 外 壁面 的 压力 系数 绝对 值 分 布 Re = 105 
Fig. 10 The wall pressure coefficient absolute value 
distribution of inner wall and outer wall, Re = 105 
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0 
s/H 
(b) 外 壁面 


图 11 外 壁面 的 压力 系数 绝对 值 分 布 Re = 108 
g. 11 he wall pressure coefficient absolute value distribution 
of inner wall and outer wall, Re = 106 


F 


ped 


< 
4 带 时 训 修正 的 清江 的 本 所 结果、 1 


CC，SKW CC 和 SSsT KW 2 安 将 用 这 6 种 


曲率 修正 后 的 灌流 模型 对 
拟 计算 , 评 佑 其 模拟 性 能 
从 图 12 可 以 看 到 ,相对 于 原 滑 流 模型 ， 带 曲率 
修正 后 的 SKE 模型 在 内 侧 壁 面 处 实验 更 加 吻合 , 而 
其 他 曲率 修正 后 的 模型 在 内 侧 壁 面 模拟 性 能 反而 变 
差 , 但 在 外 侧 壁 面 附近 速度 值 有 所 增加 ， 同 实验 值 


时 的 流动 进行 模 


吻合 较 好 ; 整体 来 看 , 对 于 速度 分 布 模拟 结果 , 曲率 。 


修正 后 的 SKE 汕 流 模型 比 原 模 型 的 模拟 性 能 略 
些 。 而 对 于 另外 5 个 灌流 模型 ， ea 时 性 

能 基本 上 没有 改善 甚至 个 别 清流 模型 还 有 一 些 恶化 
(限于 篇 幅 没 有 给 出 ) 也 就 是 说 ， i 
修正 方法 ,对 于 这 种 U 型 弯 管 流动 , 未 能 很 好 考虑 
其 中 的 湛 流 非 平 衡 输 运 机 理 和 各 向 异性 机 理 ， 没 能 
对 预测 性 能 带 来 本 质 的 改善 ， 还 需要 进一步 开展 深 
入 研究 . 

为 了 进一步 探究 曲率 修正 对 灌流 模型 模拟 性 能 

的 影响 ， 还 对 比 了 实验 测 基 截面 处 的 消 动 能 分 布 情 
况 (限于 篇 幅 未 全 部 给 出 ). 其 中 , 图 13 给 出 了 在 90" 
截面 处 SKE 模型 采用 曲率 修正 前 后 的 对 比 情况 。 修 
正 后 的 SKE 渍 流 模型 对 滑动 能 的 模拟 性 能 有 很 大 的 
改善 , 在 内 侧 壁面 处 的 庙 动 能 峰值 基本 消失 了 , 而 在 
外 侧 壁 面 处 的 峰值 有 所 增加 ， 这 与 实验 结果 趋势 相 


可， 
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吻合 . 同时 也 说 明 曲 率 是 造成 该 种 汪 动 能 峰值 出 现 
的 最 主要 的 原因 . 而 对 于 另外 5 个 滑 流 模型 , 曲率 修 
正 使 得 这 一 位 置 的 济 动 能 分 布 有 一 点 改善 , 对 SKW 
模型 影响 稍 大 一 些 ， 对 于 另外 4 个 影响 很 小 . 总体 
来 说 ， 和 预测 的 速度 分 布 类 似 ,， 曲率 修 正 对 SKE 的 
性 能 影响 最 为 明显 。 


1.0 5 
mn a EXPRe=106 
0.8 Pe SKE 
-- SKECC 
0.6 
1 
0.4 - 
i 
a 
FA 
a 
0 02 04 06 08 1.0 
u/Urr 
图 12 KE 和 SKE CC 在 s/H =2 截面 的 速度 分 布 曲 
线 ，Re = 105 


Fig. 12 The velocity distribution with SKE and SKE CC 
turbulence models, s/H = 2, Re = 106 


m mEXP Re=105 
KE 


-- SKE CC 


5 10 15 20 23 30 33 40 45 
103K/ URer 


图 13 SKE 和 SKE CC 模型 在 90° 截面 的 滑动 能 分 布 曲 
线 ，Re = 106 
Fig. 13 The velocity distribution with different turbulence 
models, 0 = 90°, Re = 106 


与 对 速度 釉 动 能 的 模拟 结果 影响 类 似 ， 曲 率 
修正 对 内 壁面 处 的 摩擦 系数 的 模拟 结果 ,也 是 SKE 
模型 的 改进 最 为 明显 (如 图 14 所 示 ， 其余 模 型 结果 
限于 篇 幅 未 给 出 ). 带 修正 后 的 模型 模拟 到 再 附 点 的 
位 置 会 比 实验 值 推 后 ,使 得 下 游 直 管 段 的 模拟 结果 
偏 小 。 对 于 内 壁面 的 静 压 力 系 数 ， 其 影响 结果 与 摩 
擦 系数 类 似 , 这 里 就 不 再 费 述 . 综合 以 上 对 比 , 表明 
SKE 模型 对 曲率 修正 最 为 敏感 . 
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图 14 SKE 和 SKE CC 模型 内 壁面 的 摩擦 系数 绝对 值 分 
布 ，Re = 106 
Fig. 14 The wall friction coefficient absolute value 


distribution of inner wall, Re = 106 
十 > 
5 结 论 


本 文 针对 二 维 U 型 弯 管 内 的 流动 , 研究 了 常用 
的 8 中 清流 模型 对 弯 管 流动 的 模拟 性 能 ， 系 统 分 析 
了 常用 涡 务 济 流 模型 对 弯 管 流动 的 模拟 能 力 ， 同 时 
用 带 曲率 修正 的 灌流 模型 对 流动 进行 模拟 , 通过 对 


ee ea dt, 


了 以 下 几 点 结论 : 


异性 较 大 .各 消 流 模型 预测 流动 从 并 科 
的 速度 比 实验 结果 慢 . 相 比 Bg 集 3b5， 各 清流 模型 
在 Re = 10" 时 对 速度 分 布 的 榴 拟 结果 与 实验 结果 更 
加 吻合 . 

2) 各 济 流 模型 对 分 离 涡 大 小 模拟 结果 偏 小 . 
SKE 没有 模拟 到 分 离 流 动 的 存在 ， 模 拟 效 果 最 差 ， 
整体 来 看 ， 当 雷诺 数 较 高 时 , RSM 模型 和 RNG KB 
模型 对 分 离 的 预测 能 力 较 好 . 


3) 从 江 动 能 分 布 和 满 流 剪 切 应 力 分 布 的 模拟 结 、 
果 来 看 ,RSM 模型 比 其 他 模型 好 。 同时 ， 对 于 摩 扼 y 


系数 和 静 压力 系数 的 模拟 结果 , RSM 模型 模拟 渴 结 
果 也 与 实验 最 为 吻合 . 

4) 曲率 修正 对 SKE 以 外 的 清流 模型 预测 性 能 
基本 上 没有 改善 ,SKE 模型 对 曲率 修正 的 影响 最 为 
敏感 。 相 对 来 说 , 修正 后 的 SKE 模型 对 于 速度 、 滑 
动能 、 内 壁面 处 的 摩擦 系数 的 模拟 结果 都 有 一 定 的 
改善 , 但 仍 有 不 小 的 偏差 

综 上 , 现 有 汕 流 模型 无 法 准确 模拟 U 型 管 流 动 ， 
采用 Spalart 和 Shur 提出 的 曲率 修正 方法 ， 也 不 能 
很 好 的 计 入 大 曲率 流动 中 的 湾流 非 平衡 输 运 特性 和 
各 向 异性 ， 针 对 大 曲率 流动 的 清流 机 理 和 汕 流 模型 
改进 工作 还 需要 进一步 开展 
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